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As it was shown earlier [Dzero, Gor'kov, and Zvezdin, J. Phys.:Condens. Matter 12, L711 (2000)] 
the properties of the first-order valence phase transition in YbInCu4 in the wide range of magnetic 
fields and temperatures are perfectly described in terms of a simple entropy transition for free Yb 
ions. Within this approach, the crystal field efi'ects have been taken into account and we show that 
the phase diagram in the B ~ T plane acquires some anisotropy with respect to the direction of an 
external magnetic field. 

PACS numbers: 75.30.Kz, 75.30. Mb, 71.23.An 



As it is well known, YbInCu4 undergoes a first-order 
valence transition at 42 K, accompanied by a small 
change in volume of the order of 0.5%. This transition 
is quite similar to the 7 — a transition in metallic Ce 
(for the phase diagram of Ce, scot). It turns out, that 
YbInCu4 is the only stoichiometric compound known in 
which an isostructurai yalence transition at ambient pres- 
sure has been observedH. However, as it was pointed out 
inp, the isostructurai valence transition is just the ex- 
treme limit of the very common competition that occurs 
between local-moment and itinerant behavior in many 
strongly correlated compounds. 

The valence transition induced by en external mag- 
netic field in YbInCu4 and its allays has been studied 
ind. The most interesting result ofEl is that the data ex- 
tracted from the resistance measurements can collapse 
all of the pressure dependent data, as well as that from 
doped variants of YbInCu4 at ambient pressure, onto a 
universal B — T phase diagram (here B is a magnetic 
field, r is a temperature). 

The physics which can be responsible for the transition 
in YbInCu4 has been discussed by several authors. One 
of the first attempts to describe the transition in metallic 
Ce, the one which is similar to the transition in YrbJnCu4, 
was the Falicov-Kimball-Ramirez (FKR) modelllJ. An- 
other approach in which the j — a transition is ascribed 
to the Mott's first-order transition ia. a subsystem of /- 
electrons has first been discussed irJiil. 

Very often the (7 — a) transition in Ce is interpret£d_iii 
terms of the Kondo Volume Collapse (KVC) modelBO. 
In the KVC model Ce atoms at the transition are treated 
as Ce'^"'"-ions in the both a and 7 phases (approximately 
one electron in the /-shell) , although in the two different 
Kondo regimes. As it is known, the Anderson impurity 
model reproduces the Kondo behavior in the regime when 
charge fluctuations are fully suppressed, and provides for 
the Tk the expression: 

Tif (xexpj-^^j, (1) 

where | e!- | is the efi^ective position of the localized level 



below the chemical potential and the levels width F c>c 
V'^i'{ep) depends on the hybridization matrix element, 
and the density of states at the Fermi level, i^{ep). The 
KVC modeltj connects the first order transition with 
strong non-linear dependence of the Kondo scale (|l|) 
(I £j |<C F) on the volume through the volume depen- 
dence of the hybridization matrix element (in Ce change 
in the unit cell volume is large, Sv/v ^ 20%!). 

Nevertheless, the KVC model seems not to be applica- 
ble in case of Y|bInCu4, where the volume changes are 
extremely smallErQ. On thatj-xeason, the FKR model 
has recently been revisited irO. It is interesting, that 
although being somewhat sensitive to the choice of the 
model parameters, the elliptic shape for the phase tran- 
sition line in the (S, r)-plane, observed inEI is preserved 
in the calculationsli3. This is probably due to the same 
mechanism as above , i.e. due to large differences in the 
energy scales for the two phases (it seems however that 
the constant a in (|^) strongly depends on the parameters 
choice) . 

As it was discussed inB, the universality of the first- 
order transition line for YbInCu4 and its alloys in the 
B ~ T plane can be described in terms of an entropy 
first-order transition between the local f-moment phase 
and another phase with a compensated moment. It 
was also suggested irJj that the mixed-valence transi- 
tion is driven by the change in the electronic screening: 
high-temperature phase can be described as a band-like 
semimetal with a small carrier concentration and accord- 
ingly screening is weak, what favors to localization of the 
/-electrons. At lower temperatures after a phase tran- 
sition occurred, even the /-electrons form a band state, 
so that a small change in occupation numbers would not 
contradict to an emergence of a large /-like Fermi sur- 
face. 

In this paper we would like to address the issue of how 
the phase diagram of YbInCu4 in the B — T plane is 
affected by taking into account the crystal field effects 
and, as a consequence, an appearance of anisotropy of 
a phase diagram with respect to the direction of an ap- 
plied field. We also would like to analyze the relation. 
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experimentally obtained i: 

ay-sBcn — 71,0 
in terms of the crystal field Hamiltonian: 

H — Hcrystal + QJl^B^ ' B. 



(2) 



(3) 

When magnetic ion is placed in a cubic environment, 
the spatial degeneracy of its angular momentum is re- 
moved by the electrostatic field due to the neighboring 
charges. For example, the J — b/2 multiplet for a Ce 
ion is split into a Fy doublet and Fg quartet while the 
J = 7/2 multiplet of an Yb ion is split into a Fg doublet, 
F7 doublet and Fs quartet. 

For the J = 7/2 the wave functions for representations 
Fg, F7 and Fg are given by ID: 




According toci the first-order transition line in the 
{B,T) plane is determined by the equation: 



T ■ S{B, T) ~ const. 



(7) 



The entropy is determined by the Yb^+ multiplet struc- 
ture only. Taking the crystal splitting effects into ac- 
count, the entropy is given by: 



T.5(i?,r) = -rin|x:exp(^-^)| 



(8) 



where An are the eigenvalues, obtained by a solution of 
an eigenvalue problem for (|^) on wavefunctions l^^- 

Now we have to find eigenvalues A„. Re- writing the 
last term in (0) as: 



gjMBJ-B = J,/?^ + J+/?_H- J_/3+, (9) 



where J± = 4 ± iJy and /3± = {B^ ± iBy), (3^ 
gjfJ-sBz, matrix elements H^j of (|^) are given by: 



V35 fj 
3 1-'^ 

^/35 o 

— —p+ 

3/3_ 

^/3+ 
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1/5- 



g/3z 

E^-lP, 3/3+ 

-^/3+ 3/3_ Es + \p.. 



(10) 



where the matrix elements of Hcrvstai are defined as: // Va \2 oqca 

I^E,-Es-l^f3A +^/3n, (12) 



< Fg|i/crysta;|r6 > — Eg, < rj\Hcrystal\^7 > — E^ , 

< rs\H,rystal\rs >- Es, (11) ^^^^(^^) ^ ^ j^E^ + Es - y/3.± 



21' / 441 



16 \^ 8960 



As it turns out, the secular equation: 

det\ \H., - XS., I H 0, \/ [^^ - ^« + i'^^ j + lif ^ 

with Hamiltonian matrix, given by (|l0|) can not be solved ^ ^ Ie E ^/3± 

exactly in its general form, but there exist the exact so- sfiKPz) 2 1 ^ ^ 7 
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lutions for particular cases, such as B = (0,0, i?z). The 

eigenvalues in that case are: / ^ ^ 16 ^ \ 768 . 

d IEt-Es + I + -j^Pt } , (14) 

Ai,2(/3.)-i|s6 + £;8 + y/3.± 
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Then 



E7 — Es 



16 



-A 



768 
l9"^ 



(15) 



In what follows we present the numerical result, when ex- 
ternal field is in the plane B = {Bx, By, 0) and analytical 
solution for the case mentioned above. 

Now, we can define a constant in Eq. (|^) as follows: 



SEe,8 



(16) 



5, 



In 



4 + 2 exp 



According toEl SEq,s 
Ej- Eg- 3.8meV. 



Er — Ea 



2 exp 



5Ej^ 



3.2meV and (5£'7,. 
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FIG. 1. The phase diagram for the YbInCu4 shows some 
anisotropy with respect to the direction of an external mag- 
netic field: B is along one of the main cubic axes (dashed 
line), B is in the easy plane (solid black line) and B is along 
the one of the main cubic diagonals (solid gray line). 



As we already mentioned, the equation (16) defines a 
phase diagram in the B — T plane. The results of our 
calculation are plotted on Fig. 1. As we can see, there 
is a strong anisotropy in a phase diagram with respect 
to the direction of an external magnetic field. We also 
have calculated magnetization as a function of an exter- 
nal magnetic field for a given temperature. Fig. 2. As it 
turns out, magnetization also depends on the direction 
of an applied field. 

In the rest of the paper, we would like to write down 
an equation for the phase boundary, when B = (0, 0, Bz). 
Let us introduce the following notations: 



b = B/B,n, T = T/r,o, tan((^) =b/T {0<^< -). 



b = U^'^if) tan(i^), 
1 



^(o,7;o) 



^exp -A„(v3) 



(17) 
(18) 

(19) 



with Xn{f) being the eigenvalues ( |12[ - |15| ). Thus, the 
equation for the phase transition line is given by: 



P+T^ = i?((^), R{cp) = cos{ip)U{ip) 



As we see on Fig. 3, deviations of R{(p) from 1 do not 
exceed 10%. Generally speaking, as one may see from 
Table 1, the a value is relatively robust with respect to 
including new type of interactions (or anisotropy) in our 
model. For example, if one takes the susceptibility of a 
lower phase into account, it decreases a by reducing the 
value for the net moment /i which in our case is less then 
AfiB and the latter corresponds to the free ion model, but 
on the other hand it increases one by reducing a change 
in the entropy at T = T^q. As we see from results of our 
calculations of Bco (Table 1), when only crystal fields 
effect are taken, the change in a is less then 10%. 
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FIG. 2. Magnetization curve shows some anisotropy with 
respect to the direction of an external magnetic field in agree- 
ment with our result for the phase diagram where the same 
type of anisotropy is found. Calculations were performed for 
T = 0.75T„o 

We have to mention that although the crystal field ef- 
fects do not change a value of a parameter in (0), it is not 
the case if one will try to verify the relation (^) in terras 
of exact solution for the Kondo model with J = 7/20. 
As it turns out, in that case, the value for a strongly de- 
pends on the energy scale, given by Xfc. This result is in 
fact similar to the result obtained byllJ using the Dynam- 
ical Mean Field Theory (DMFT) approach based on the 
FRK model. Thus, experimentally verified exsistancc of 



3 



anisotropy in the phase diagram will serve as a proof of 
our initial idea of free Yb ions model. 




FIG. 3. Deviations of R{ip) from 1 do not exceed 10%. 

To summarize, we have shown that the phase diagram 
for the first order valence transition in YbInCu4 in the 
B — T plane acquires some anisotropy with respect to 
the direction of an external magnetic field if crystal-field 
split multiplets are taken into consideration. We have 
also found that within the present framework, that the 
anisotropy of the critical field value Bco is of the order of 
~ 2 Tesla and in principle can be seen experimentally. 

TABLE I. Values of Bco and a obtained from the phase 
diagram Fig. 1 for different magnetic field orientations. 



Field Orientation 


Bco, Tesla 


a 


[001] 


30.9 


2.29 


[111] 


30.1 


2.34 


[110] 


31.95 


2.24 



Preliminary experiments have been carried out in the 
National High Magnetic Field Laboratory (Tallahassee) 
and showed relatively good agreement jssdth theoretical 
prediction in terms of the present modelEj. 
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bringing this problem to my attention and for very use- 
ful comments on this paper. I also thank P. Schlottmann 
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exact numerical solutions for the asymmetric Anderson 
model and S. Nakatsuji for useful discussions. This work 
was supported by the NHMFL through the NSF cooper- 
ative agreement DMR-9527035 and the State of Florida. 
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